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Abstract 

Back-to-Back Correlations of particle-antiparticle pairs are related to the in-medium mass- 
modification and squeezing of the quanta involved. They are predicted to appear when hot and 
dense hadronic matter is formed in high energy nucleus-nucleus collisions. The survival and mag¬ 
nitude of the Back-to-Back Correlations of boson-antiboson pairs generated by in-medium mass 
modifications are studied here in the case of a thermalized, finite-sized, spherically symmetric ex¬ 
panding medium. We show that the BBC signal indeed survives the finite-time emission, as well 
as the expansion and flow effects, with sufficient intensity to be observed at RHIC. 
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I. INTRODUCTION 


Recently, it has been shown [1, |2] that large back-to-back correlations (BBC) of particle- 
antiparticle pairs of bosonic particles might appear in high energy nucleus-nucleus collisions 
as a consequence of in-medium mass modification of the bosons. Detailed calculations in¬ 
dicate that the BBC signal app ears for values of transverse momenta below 1-2 GeV/c. 
More recently, it was shown [3] that BBC of similar strength might appear for fermionic 
particles as well. The main physical ingredient used in the evaluation of the effects of 
in-medium modified masses on two particle correlation functions is a quantum-mechanical 
correlation induced by a nonzero overlap between in-medium and free states. The induced 
quantum mechanical correlation can be represented in terms of two-mode squeezed states of 
the asymptotic, observable states and is implemented through a Bogoliubov-Valatin trans¬ 
formation. 

The possibility of measuring a significant BBC signal in heavy-ion collisions opens new 
interesting possibilities for accessing the properties of the matter formed in such collisions. 
The BBC signal is linked to in-medium mass modifications of hadrons in the hot and dense 
environment the detected particles experience before freezing out and in this sense BBC mea¬ 
surements provide independent pieces of information on medium modifications from the ones 
obtained from dilepton yields and spectra. However, there are several additional physical 
effects that interfere with mass modifications of the detected particles in the interpretation 
of the BBC signal. 


. [H, 0] and 


All studies in Refs. Ill, \M and |3|| were restricted to infinite, static media. For testing the 
robustness of the effect, we generalize the previous studies to a more realistic situation of 
mass modification in a finite-sized, expanding thermalized medium. In this first, investigation 
of such effects, we use a simple hydrodynamical model for the expansion and simple three- 
dimensional Gaussian profile for the size of the system. Although simple, the model is rich 
enough to indicate the influence of the expansion and finite-size effect on BBC. 

In the next section we review the basic ingredients of the model and in Section IITII 
we generalize the model to include expansion effects. Numerical results are presented in 
Section Hv| and Section m presents our Conclusions and Future Perspectives. 
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II. REVIEW OF THE MODEL - INFINITE HOMOGENEOUS MEDIUM 


In the present paper we concentrate on bosonic BBC and restrict the discussion to cases 
where the boson is its own antiparticle - like the (j) meson. We are interested in the two- 
particle correlation function 


C 2 (k!,k 2 ) 


W( k^A^k,)’ 


( 1 ) 


where iV^kj) and iV 2 (k 1 ,k 2 ) are, respectively, the invariant single-particle and two-particle 
momentum distributions 


d 3 N 


iVi(ki) = = a; kl (a^a^), 

V 2 (ki,k 2 ) 


dkx 

(UkiCUkj (®ki®k2 < ^k2®ki) 

= Co , k 1 ^k2 ( a ki®ki) (0'k 2 a k 2 ) “I - (^kj 0 ^) ( a k 2 a ki) 


( a ki a k 2 ) ( a k 2 a ki) 


( 2 ) 

( 3 ) 


where a^. and ak are free-particle creation and annihilation operators of scalar quanta, and 
the angular brackets mean thermal averages. The factorization of the expectation value 
of four operators into products of expectation values of two operators in Eq. © has been 


derivec 
Refs, y 


asji generalization of Wick’s theorem for locally equilibrated (chaotic) systems in 
. Introducing the chaotic and squeezed amplitudes as 


^c(ki,k 2 ) — A/Md^k2 ( a ki a k 2 )> 
^(ki,^) = ^/cUki^ka (Q'ki®k 2 )5 


( 4 ) 

( 5 ) 


the two-particle correlation function can be written as 

|G c (k!,k 2 )| 2 


C 2 (ki, k 2 ) — 1 + 


+ 


|G s (k 1 ,k 2 )p 


( 6 ) 


G c { k l5 k 1 )G c (k 2 , k 2 ) G c ( k 1} k 1 )U c (k 2 , k 2 )‘ 

The G c (l,2) is the usual Hanbury-Brown-Twiss (HBT) amplitude and G s (l,2) is BBC 
amplitude. 

The thermal average of an operator O, ( O) = Tr(pO), is calculated with a density ma¬ 
trix p corresponding to medium-modified, thermalized quanta. The crucial point is that the 
in-medium thermalized quanta are not the ones detected. The detected quanta have energy- 
momentum k 11 = (cuk, k), a>k = k 2 + m 2 and are described by the creation and annihilation 
operators al and a k- However, if we denote by b £ and b k the creation and annihilation oper¬ 
ators of in-medium, thermalized quanta with k*^ = (fik; k), = k 2 + m 2 (|k|), we can relate 
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the (a^, (ik) to (bl,bk) through a Bogoliubov-Valatin (BV) transformation. Specifically, the 
annihilation operator a k x for the asymptotic quanta with momentum ki is related to the 
in-medium operators b kl and as |l|: 

^ki = Ck x &ki + ■‘>h k , &^ki = Ci + S-i, (7) 

where we have introduced the notation C\ = bk L and S -1 = s* k ^b k| to simplify later 
notation, and 

c k = cosh[/ k ], s k = sinh[/ k ], fk = \ lo S • ( 8 ) 

The BV transformation for the creation operator a kl is obtained from Eq. 0 by Hermitean 
conjugation. As is well known, the Bogoliubov transformation is equivalent to a squeezing 
operation, and this motivates calling fk the mode-dependent squeezing parameter. In this 
way, it is the squeezing parameter /k that carries the in-medium effects. Using the BV 
relation, we obtain for the thermal averages in Eqs. 0 and © 


C c (k 1 , k 2 ) = [(C\C 2 ) + (5'_ 1 5'1 2 )] , 

C s (k l5 k 2 ) = [(^ 2 ) + (CiSU) ■ 


(9) 

( 10 ) 


After performing the thermal averages indicated above, with the help of a thermal den¬ 
sity matrix p corresponding to the in-medium modified, thermalized quanta, the resulting 
expressions for the case of an homogeneous medium are 


C c (l,2) = | 
G s { 1,2) = { 


Ci, 2 
(2t r) 3 
Ci,2 
(2tt) 3 


|ci, 2 | 2 ni, 2 + |s_i,_ 2 | 2 (n_i _2 + 1)J | V<$i )2 , 

s* 1,2 C2-1 n_i,2 + Cl _2 S *_ 2)1 (ni -2 + 1) } V5 \- 2 , 


( 11 ) 

( 12 ) 


From Eq. HU and (H2D it is easily seen that, in the approximation of a sudden freeze 
out, and in the case of a homogeneous medium, G c (ki, k 2 ) oc V<5i, 2 and G s (ki, k 2 ) oc V8\ _ 2 . 
Therefore, the amplitudes C c (ki,k 2 ) and G s (ki,k 2 ) are non-vanishing only for ki = k 2 and 
ki = —k 2 respectively. I 11 the expression for the two-particle correlation function the volume 
factors cancel out, and we obtain jfj 


C 2 (k,k) 


C 2 (k,-k) 


2 , 

. |c k gk n k + C- k 5h k ( n ~k + 1) | 2 
ni(k)ni(-k) 


(13) 

(14) 
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where ni(k) is defined by 


(15) 


with 


Afi(k) 


V 

(2t r) 3 


07 k ni(k), 


ni(k) 


Ck| T |s_ k | (rm k H - !) , 


(16) 


and n k is the Bose-Einstein distribution function of the in-medium quanta with energy 
f2 k at temperature T. The exact value of the intercept, C^k, k) = 2, is a characteristic 
signature of a chaotic Bose gas without dynamical 2-body correlations outside the domain 
of Bose-Einstein condensation. 

We should note that Eq. m is valid only in the rest frame of the medium, i.e., the 
correlation is back-to-back only in the rest frame of the matter. In the next section we extend 
the model to a medium with finite size corresponding to a fireball, which is exploding with 
a position dependent flow velocity field distribution, so that only the central point of this 
exploding fireball is at rest in the frame of the observation. 


III. SPECTRA AND CORRELATIONS FOR MASS-SHIFTED BOSONS IN FI¬ 
NITE EXPANDING SYSTEMS 


We are mainly interested here in the study of the squeezed correlation function - first and 
third terms of Eq. ©• For studying the expansion of the system we adopt for the emission 
function the non-relativistic hydrodynamical parameterization of Ref. jj, which was shown 
later to actually be a non-relativistic hydrodynamical solution. In this model the fireball 
expands in a spherically symmetric manner with non-relativistic four-velocity = 7 ( 1 , v ), 
with 7 = (1 — v 2 ) -1 / 2 « 1 + v 2 / 2 , where 

v = (u)r/R, 

(u) and R are, respectively, the mean expansion velocity and the radius of the fireball. Thus, 
we divide the inhomogeneous medium into independent cells and assume that Eqs. © and 
m can be evaluated locally within each cell using the BV transformation of Eq. © - and 
its Hermitian conjugate. Then, the amplitudes G c and G s can be written in the special 
form derived by Makhlin and Sinyukov 4J, which are given by Eqs. (22) and (23) of Ref. Q], 
namely 


G ? c(kl,k 2 ) — J /27r) 3 ^ e * 91 ’ 2 X { I C l,21 2 '^1,2 (^) + 1^—1,—21 2 ['^-1,-2 (x) + 1]}; 


(27r 


(17) 
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G s (ki,k 2 ) = 


/ ^^’ 2e2lKl 2X { S *~ 1,2C2 '- in ~ l ' 2 ^ X ) +Cl ~ 2 ' s - 2 - 1 i ni - 2 ( x ) + 1 ]}- 


( 18 ) 


Here dra^ix) = g? 3 E M (x; r/) F(T/)dr/ is the product of the normal-oriented volume element 
depending parametrically on the freeze-out hypersurface parameter 77 and on its invariant 
distribution function A(t/). We should notice that, in the particular case in which each 
dra^x) of Eqs. (17) and (18) is parallel to it 11 , that is, the emission from an elementary cell 
mentioned above occurs instantaneously in its proper frame, the exponential factor there will 
give rise, upon integration over the cell assuming it is large enough, to the same factor 
or 81-2 that were present in Eqs. (11) or (12), respectively. As mentioned at the end of Sec. 
II, the arguments of 5 here are not ki and k 2 of the left-hand side, but should be understood 
as given in the proper frame of the cell. In what follows, the condition of instantaneous 
emission in the proper frame of each cell is assumed to be approximately verified, since our 
calculation is non-relativistic. However, we should remark that due to the fact that our 
elementary cells are not always large, the correlation described above is only approximately 
back-to-back. 

The other quantities appearing in Eq. m and m are n ly j(x) = n(x,Kij), the local 
density distribution, and c hJ = cosh[/jj(x)] and s hJ = sinh[/jj(a;)], squeezed functions with 


fi,A x ) = 2 log 


AT,-(a;) u^x) 




K*J(x) u v {x) 


(19) 


where u^^x) is the local flow vector at freeze-out. The relative and the average pair four- 
momentum coordinates are defined as g° 2 = uq — cu 2 , qi )2 = ki — k 2 , AT ° 2 = (uq +u; 2 )/2, and 
K 1)2 = (ki + k 2 )/2. Also, we identify in-medium and squeezed quantities by superscripted 
asterisks. The relative ( x ) and total four momenta A' M of particles 1 and 2 are given by 


1 ij(x) = k?(x) - kf(x), 




X) 


k?(x) + k$(x) 


( 20 ) 


where kf(x) for i = ±1, ±2 are given by 


= Xi kl (x) u"(x) ± kf(z), w kl (l) = y m 2 - k, ki„ = kfu„(x). 


( 21 ) 


with kf orthogonal to u^(x)\ 


ki = k!^ — ki ■ u(x) u fl (x ), 


( 22 ) 
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The corresponding in-medium quantities are given by 


1 r 


lij(x) = k‘"(x) - k]“(x), Kij(x) = - [kT(x) + k‘k(x) 


(23) 


and 


k%(x) = Sl* l (x)u“(x)±k?‘(x), 


k) - krk„ = OvM, (24) 


with 

= k* 11 — k* ■ u(x ) u^(x). 


(25) 


Now, it is not difficult to show that k*^(x) = k^(x) and therefore no star is necessary in k for 
the in-medium quantities. It should be noted that this equality was not clearly emphasized 


f 0. 


in Ref. 2||. An important aspect of these relations is that due to the mass modification, the 
energy in the local co-moving frame is modified from uJ^Jx) = fcf (x) u^(x) = k± t (x) (x) 
to = K^{ x ) u ^ x ) = kTi{ x ) u v( x )j without modifying the component of the four- 

momentum orthogonal to the four-velocity. The above definitions are the detailed write- 


0 > 


up of similar definitions of Ref. [2||, where a more succinct notation has been used and 
the misprint signs (=F, instead of ±) in Eq. (27) and (28) of Ref.|2| have been corrected, 
respectively, in expressions for k^^x) and k^^x) above. These definitions of momenta are 
illustrated in Fig.|2 corresponding to the relativistic and non-relativistic limits, in parts (a) 
and (b), respectively. 

Using the above expressions, the squeezing parameter can be evaluated as 




lj \ K ij( x ) M x ) 
2 ° to _K%{x) u v (x) 


1 j r ^Qj) +u^(x) 

2 ° g _U kl (ic) + ^kj(^) 


f±i,±j{ x )i 


(26) 


where the spatial coordinate dependence enters through the position dependence of either 
the four-velocity or the in-medium mass modification, or both. However, it does not matter 
which of the locally back-to-back momenta are used for the evaluation of the amount of 
squeezing. 

In this paper, we focus on the effects of expansion on the back-to-back correlations - or, 
with other words, does the flow wash out the signal for these correlations or not? Although 
the formalism of the squeezed back-to-back correlations was worked out with expanding 
systems in Ref. 2, no detailed investigations were performed to quantitatively study e.g. the 
strength of the signal with varying the strength of the flow. 
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FIG. 1: The definition of the various momentum possibilities discussed above (where the index 1 
was chosen for simplicity), and the illustration of the dependence of the back-to-back momentum 
pairs on the four-velocity and on the mass-shift, is shown in (a) using relativistic notation. Solid 
lines represent the locally back-to-back momentum pair and its components for asymptotic quanta; 
dashed lines represent the same for in-medium modified, mass-shifted quanta. In (b), the analogous 
situation is illustrated using the non-relativistic notation. 


Here we intend to investigate this question in one of the simplest geometrical cases. For 
the sake of clarity, we evaluate the flow effects for a non-relativistically expanding, spherically 
symmetric fireball, that freezes out at a constant temperature T and has a Gaussian density 


profile. In this sense, we adopt the model emission function of Ref. |, that was developed to 
study single particle spectra and Bose-Einstein (HBT) correlation functions in the simplest 
possible case of expanding systems. Later on it has been realized, that this emission function 
corresponds to the simplest member of a new family of exact solutions of non-relativistic 
hydrodynamics jsj], which can be generalized in a straightforward manner to cylindrically 


Q, 


and cllipsoidally symmetrical expansions, as well as to the case of relativistic expansions 


and in all cases, systems that expand with inhomogeneous temperature profiles &Q- 
Before investigating in detail the effects of various kind of inhomogeneities in the flow 
profiles and in the temperature profiles, let us turn our attention to the non-relativistic limit 
adopting the simplest possible scenario for the expansion that leads to analytic forms [7, Q| . 

For the sake of clarity, we present the explicit expressions in the non-relativistic limit 
of the above quantities, which is the appropriate limit for our non-relativistic flow model. 
Writing u M = 7 ( 1 , v) and using 7 = (1 — v 2 ) -1 / 2 ps 1 + v 2 /2 « 1 , we have that Eq. 
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leads to 


and therefore 


k i — mv = k* — m*v, 


(27) 


k±i(r) 

= m v(r) =)= m v(r) ± k, : , 

(28) 

k ii( r ) 

= m* v(r) =F m* v(r) ± k*. 

(29) 


With this, we obtain 

K‘j(x) u F ss m, + T- |[Ky - m,v(r )] 2 + i(q*,) 2 ), 


(30) 


where the total and relative local in-medium momenta are given by 


Kt i2 ( r ) = + k 2 ) + i 

(m* - m) v(r), 

qt 2 ( r ) = 

: (k X - k 2 ) 

(31) 

K i,-2( r ) = m * v (. r ) 

+ 2^! _ ^ 2 ) ’ 

qi,- 2 ( r ) : 

= -2mv(r) + (ki + k 2 ) 

(32) 

K* i >2 (r) = m*v(r) 

- |(k L - k 2 ), 

q-i, 2 ( r ) : 

= 2mv(r) - (ki + k 2 ) 

(33) 

-i,_ 2 ( r ) — ( <m * + m)v( r) 

— y* 21+ ’ 

q-i,- 2 ( r > 

'k 

1 

K 

to 

(34) 


The unstarred Ky and qy momenta are obtained from the above by replacing m* by m in 

these expressions. Note that q* • = q y, as it should be. These relations imply, that in the 
non-relativistic limit, 

Kb — m*v = Ky — mv. ( 35 ) 

In discussing finite-size effects, we distinguish between the volume of the entire thermal- 
ized medium, denoted by V, and the volume filled with mass-shifted quanta, denoted by V s . 
Naturally, V s < V in the general case. In the derivation of the expressions for G c ( 1,2) and 
G s ( 1,2), for simplicity, we introduce a Gaussian profile function in the integrands, i.e., we 
consider that the volumetric region where the mass m* is significantly modified is smooth 
and Gaussian in shape. In other words, instead of considering a particular domain of inte¬ 
gration, we perform the spatial integrals for G c (l,2) and G c (l,2) using a Gaussian weight 
e ~r-/ 2 Rl j n ^ ie integrands, extending the integration region to infinity. Specifically, we have 
for G c (l,2) and G s (l,2) 
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( 36 ) 


G c (ki,k 2 ) 


G a (ki,k 2 ) 


El,2 

Wi 


j d 3 re- i(kl - k2) r | e - r2/2i? ' ^|c(l, 2 )\ 2 n\^{ x ) 


+ |s(—1, —2)| 2 [n* 1 _ 2 (a:) + l] j + (l - e r2/2jR ") ni )2 (a;)|, 

J d 3 r e -*(ki+k 2 )T e -r»/ 2 ^ 2) c(2, -1) n* li2 (x) 


(2vr) 3 

+ c(l, —2) s*(—2,1) l + <_ 2 (x) 


(37) 


where n*j(x) means that the local distribution function is to be evaluated with in-medium 
momenta, i.e. n*j(x ) = n(x, K*- ). The integral over the factor (l — e - r2 / 2i ^ represents the 
integration over the region where there is no mass shift, corresponding to the region V — V s . 
In this region, we have that the squeezing factors become c(i,j ) —> 1 and s(i,j) —> 0. 

In order to proceed, we need the expressions for n*(i,j) and n(i,j). We consider their 
Boltzmann limit, 


n(x, K t .j) ~ exp 


K ti u A x ) - aM 

T(x) 


(38) 


and the same for n*j(x ) with K h] (x) replaced by K*j(x). Considering that the chemical 
potential in the model of Ref. can be written as fi(x)/T(x ) = n 0 /T — r 2 /2R 2 , and making 
use of Eq. (13011 . it is easy to show that 


n 


1,2 


(r) = Uq exp 


[(ki + k 2 )/2 — m(u)r/R\ 2 (k x - k 2 ) s 


2 R 2 2m*T 

w - 1 ,- 2 ( r ) = n -i,2( r ) = nl r2 ( r ), 


8 m* T 


(39) 


where 



(40) 


This factor is proportional to the mean multiplicity, and can be determined in principle 
from the absolute normalization of the single particle spectra. The corresponding unstarred 
nij(r) are obtained from n*j(r) by replacing m* by m in Eq. m- 

When evaluating the spectra and the correlations from this model, we realize that a 
mathematically equivalent problem has already been considered in Ref. [7]. By replacing 
m —> m* and t 0 — > Rqtti^/( u)m in the equations of Ref. Q, the results obtained there can 
be directly transcribed here. 

Due to the equality in Eq. (1331) . we see that the accounting for the squeezing effects can 
be simplified for small mass shifts (m* — m)/m <C 1, such that the squeezing parameter can 
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be written as 


= ^ lo g 




_K%(i,j,x)u v (x)_ "> log ^~~^* ^ 

The neglected terms are order of (kinetic energy/mass) 2 (masshift/mass) 2 and hence are of 
fourth order in small quantities. This limit is important, because in this case the coordinate 
dependence enters the squeezing parameter / only through the possible position dependence 
of the mass-shift and so the flow effects on the squeezing parameter can be neglected. In 
principle, the position dependence of the mass shift can be calculated from thermal held 
models in the local density approximation. Therefore, in an approximation that the position 
dependence of the in-medium mass is neglected, the c(i,j) = Co and s(i,j ) = so factors can 
be removed from the integrands and all what remains to be done are Fourier transforms of 
Gaussian functions. The final result for G r and G „ can be written as 


G c (k!,k 2 ) = 


-^ 1,2 

w> 


n* 0 (|c 0 1 2 + |s 0 | 2 ) -^ 1 , 2 (S s , R, T, m*) + |s 0 | J /p 2 (-R s , oo, oo, m*) 


^(ki, k 2 ) = 


+ n 0 Ii 2 ( oo, i?, T, m) — n 0 I^ 2 (R s , R, T, m) 
Sy 2 


(2vr) : 


c 0 | s 0 | 2 nl /i i2 (-R s , R, T, m*) + /p 2 (i? s , oo, oo, m* 


where 2 and /* 2 are the resulting Fourier integrals 


Ii, 2 (R s ,R,T,m w ) = (27T/9 2 *)) 3 / 2 exp-! 


(k? + k 2 


4 m(*)T 




(ki - k 2 ) + * 


,m(u)(ki 


2 mt*\TR 


Il 2 (R s ,R,T,m *) = (2vrp 2 ) 3 / 2 


exp < 


h*) 
(k 2 + k 2 ) 
4 m*T 


p: 

2 


1 + i 


. m(w) 


2mJT R 


(42) 

(43) 


(44) 


(k, + k 2 ) >,(45) 


with 


P(*) 


1 

S 2 


T( 1 + ^' 

S V m(*)T / 


(46) 


Finally, we include finite-time emission effects in a schematic way using for the invariant 
function F(t), that appears in the expression for d 4 cr^(x), the following expression 


F(t) = ^ T T °- e -( r “ r o)/ At ; 


(47) 


where At is a free parameter. This hnishes the derivation of all the ingredients needed to 
evaluate C 2 (ki,k 2 ), Eq. (JOJ). In the next section we will present numerical results. 


11 





















IV. NUMERICAL RESULTS 

We present numerical results for two situations regarding the volumes over which mass 
modification occurs. In the first situation the mass shift occurs over the entire volume of the 
expanding system, i.e. V s = V. This amounts to removing the factor e - r2 / 2 #! in Eqs. (ESI) 
and (El or, equivalently, take R s —> oo in I c and I s . In the second situation, we consider 
V s < V. In order to comply with the non-relativistic nature of the expansion model used 
in this paper, we present numerical results for the Back-to-Back Correlations of a 0 meson 
pair. In free space, the </> meson mass is = 1020 MeV. 


A. Bulk decay of a volume filled with mass-shifted quanta 


In this case, we suppose that the mass-shift occurs in the entire volume of the system, 
for simplicity considered here as a Gaussian with r.m.s. radius R. We will focus on the 
BBC correlation function, whose generic expression consists of the first and third terms 
on the right-hand-side of Eq. ED- The detailed expressions for the amplitudes are given 
in Appendix El In what follows, we will concentrate on the value of momenta of the 
participant pair that maximizes the BBC signal, i.e., the case in which ki = — k 2 = k. The 
BBC correlation function can then be written as 


C^(k 1 ,k 2 )=CflB 0 (k,-k) 


1 | |Gk(k, —k)| 2 
G* c (k, k)G c (—k, —k) 


\GsQt, -k)| 2 

[Gc( k, k)] 2 


(48) 


In the above equation, we have used the fact that the single-inclusive distribution, G c (k, k) = 
G c (— k, —k) depends only on the absolute value of the momentum, as can be seen in Eq. (1A8|) . 
With the aid of this equation, as well as of Eq. M) , the expression of the BBC correlation 
function is finally written as 


Cbbc( k, —k) — 1 + |c 0 sol 


2n* 0 Rl exp(—k 2 /2m*T) + R 3 


n* Q Rl (|c 0 1 2 + |s 0 | 2 ) exp(—k 2 /2m*T*) + -R 3 |s 0 | 2 


(49) 


where 


1 1 / m 2 (u) 2 ' 


R: R 2 


mJT 


m m m , , 9 

T* = T H- (u) , 

m* 


(50) 


as given in Table |U of Appendix El The quantities R t and T* are, respectively, the flow- 
modified radius and the flow-modified temperature of the system where the mass-shift occurs 
in the entire volume. 
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FIG. 2: The Back-to-Back Correlation is shown as a function of the shifted mass m* on the left 
panels, and as a function of both m* and the momentum of each particle (ki = —k 2 = k) on the 
right ones, for R = 7 fm/c, T = 140 MeV and At = 2 fm/c. The plots (a) and (b) illustrate 
better the behavior of the BBC signal seen on the right panel, when |k| = 500 MeV/c and for 
|k| = 1000 MeV/c, respectively. In both cases, the dashed curve corresponds to (u) = 0 and the 
solid curve, to ( u} = 0.5. In (c), it is shown the 3-D plot corresponding to the no flow case, with 
( u ) = 0, whereas in (d), a radial flow with v = (u)r/R = 0.5 was considered. 

In Fig.|2]we present the results for C BBC for T = 140 MeV, R = 7 fm and a finite emission 
time At = 2 fm/c, for two different flow values. We clearly see in this figure that the presence 
of a radial flow causes the BBC correlation to be higher than the case with no flow in the 
low momentum region and for the same values of m* and ki = —k2 = k, but it grows more 
slowly than in the no-flow case for increasing values of k and same m*. We also see that for 
a flow of ( u) = 0.5, the BBC signal increases for values of the momenta |k| <■ 1000 MeV/c, 
but the no-flow case surpasses the previous case for |k| > 1000 MeV/c. This conclusion is 
more easily achieved by looking into the right panel of Fig. [21 The inversion of the BBC 
behavior for that value of k roughly coincides with the limit of applicability of our non- 
relativistic approximation. This result is very stimulating, suggesting that we would have 
bigger chances of observing the BBC effect in the lower |k| side of the the BBC x m* x |k| 
region. 

In Fig. 01 we present results for different combinations of temperatures and flow velocities 
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FIG. 3: The effect of a finite emission interval on the Back to Back correlation function, as 
compared to instant emission, is illustrated in the left panel. The dashed curves have been reduced 
by a factor of 400, and the solid curves correspond to the suppression by a finite emission duration, 
of about At ~ 2 fm/c. The plot in (a) shows this effect in the absence of flow, already discussed 
in our previous paper. The plot in (b) shows the corresponding result when flow is included, with 
<u>= 0.5 (the other parameters adopted to produce the curves are R = 7 fm/c, T = 140 MeV and 
At = 2 fm/c. In part (c), we fixed the other parameters and study the influence of increasing values 
of |k| (the dashed curve corresponding to |k| = 500 MeV/c and the solid one, to |k| = 500 MeV/c). 
Finally, the plot in (d) considers the variation of the BBC curve for increasing values of |k|, as 
indicated by the arrows, keeping all the other variables fixed to their values above. 

and emission times. We see again that flow enhances the BBC correlation function for small 
values of |k|. The effect of the temperature is such that the BBC signal is stronger for lower 
values of T, and decreases the signal for higher values of T. 


B. Decay of a volume partly filled with mass-shifted quanta 

Here we consider the scenario in which the mass shift of the bosons occurs in part of the 
system volume only. In this case, the expression for the BBC correlation function is more 
complex, even in the simpler limit of the maximal effect, i.e., when ki = — k 2 = k. The 
expression of Cg can again be obtained from Eq. (14811 , but this time we must replace the 


14 







amplitudes by the expressions of Eqs. (IB 101) and (IB111) of Appendix |B] in the limit that the 
particles are back-to-back. In this case, we must use for the squeezing parameters in the 
region where there is no medium modification their appropriate limiting values 


lim 

ra *—>m 



With this, we obtain after a long but straightforward calculation the expression 


(51) 


C bbc (k, k) — 1 + 


Co So 


2 n * 0 pi expl 


^r Rs 


\n * 0 pi (|c 0 | 2 + |s 0 | 2 ) exp 


x 

2\l,2 


(R 2 + R 2 ) k 2 , ,, l9 

1 ' +i£a 0 2 


+ no TZ 3 exp 


k 2 \ 


2 mT 


J 


2m* (R 2 T + R%) 

/ ( R 2 + R 2 ) k 2 


nop exp 


2m{R 2 T + R S 2 T) 


(52) 


where the parameters T = T + m(u) 2 , 7W 2 = R~ 2 (1 + m(u) 2 /T), p ~ 2 = 1Z ~ 2 + i?“ 2 , 
p ~ 2 = R~ 2 (1 + m 2 (u) 2 /m*T) + R~ 2 are given in Table HH of Appendix |B| The two are, 
respectively, the flow-modified temperature and the flow-modified radius in the region of no 
mass-shift. On the other hand, p and p * are effective radius parameters corresponding to 
the no mass-shift region and the region where the mass-shift occurred, respectively. We see 
that they are functions of the flow parameter ( u) and the parameter R s , which corresponds 
to the radius of the mass-shift region. The parameter T* = T + m 2 (u ) 2 /m*, is the same as 
before, also written in Table [Q of Appendix f~YI 

Similarly to Eq. (©, the above Eq. m has also very interesting limiting cases. The 
first one is the case of vanishing squeezing, m* —> m, which implies that | s 0 —» 0 and 
| Co | —> 1, and the squeezed Back-to-Back Correlations vanish. 

The large momentum limit is also very interesting. In this case, the exponential, thermal 
contributions disappear, and the surviving terms come from the decay of the squeezed vac¬ 
uum to the asymptotic quanta. Both the numerator and the denominator of Eq. (E2) will 
be proportional to the square of the squeezed volume, hence Eq. (15211 and Eq. (JUJ) will be 
reduced to a form similar to 

C(k, —k ) = 1 + | Co/ Sq | 2 


which has no upper limit, and diverges for small but non-vanishing amount of squeezing, 
where |cq| —*► 1 and |s Q | —> 0. This property, the unlimited strength of the squeezed BBC-s 
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FIG. 4: The plots in this panel are similar to the ones on Fig. |2J The main difference is that here 
we assumed that the mass-shift occurred only in a smaller part of the system volume. Here, the 
back-to-back correlation is shown as a function of the shifted mass m* on the right panel, and as 
function of both m* and the momentum of each particle (ki = —k2 = k), on the left ones, for 
R = 7 fm/c, R s = 5 fm/c, T = 140 MeV and At = 2 fm/c. The plots in parts (a) and (b) illustrate 
better the behavior of the BBC signal seen in parts (c) and (d), for |k| = 800 MeV/c and for 
jkj = 1000 MeV/c. In both cases, the dashed curve corresponds to (u) = 0 and the solid curve, 
to (u) =0.5. In the 3-D plot in (c), no flow ((u) = 0) was considered, whereas a radial flow with 
v = (u)r/R = 0.5 was included in the (d) plot. 

- even if the mass modification does not happen in the whole volume - makes it worthwhile 
to look for these effects experimentally as signals of in-medium mass modifications. Again, 
for large in medium mass modifications and large momenta, the strength of the squeezed 
BBC-s will be similar to that of the HBT effect: 

C{k,-k ) —> 2 

if |c 0 |, |s 0 1 —»> oo, as in this limit, |c 0 |/|s 0 | —»► 1. 

The single particle spectra also behaves interestingly in these limits, which is discussed 
in Appendices El and |B] 

In Fig. Q] we show the BBC correlation corresponding to the hypothesis that the mass- 
shift occurred only in a smaller part of the system volume. We see a very close similarity 
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to the results corresponding to the mass-shift occurring in the entire system volume, shown 
in Fig. [21 The major difference between the two of them is that in Fig. 0] the correlation 
signal is lower, as expected, since a the mass-shift occurred in a smaller volume in this case. 


V. CONCLUSIONS AND FUTURE PERSPECTIVES 

In this paper we discussed the effects of the system expansion and flow on the back-to- 
back correlation, also limiting the system to a more realistic finite size. For simplicity, we 
restricted our analysis to the non-relativistic domain. In our study, we have also considered 
that the flow effects on the squeezing parameter were negligible. For simplicity, we have also 
assumed a 3-D Gaussian profile for the system. We showed here the effects of the decoupling 
temperature on the BBC signal, fixing all other parameters, as in the bottom right plot of 
Fig. H In this case, we observed that the BBC signal survives stronger if the decoupling 
temperature is lower. Fixing T and the other parameters, we also showed the effect of 
increasing momentum on the signal survival. More striking, we showed in Figures |21 and 0 
the conclusion about the best region for looking into the BBC effect in the m* x |k| plane: 
the search for the signal is more pronounced in the small |k| region, when we are to take into 
account the system expansion and the presence of moderate to strong flow. For higher values 
of |k|, the BBC signal would be more pronounced if the system flow could be neglected. In 
any case, what remains as the most encouraging point coming out of our present study is 
that the BBC seems to survive with measurable intensity in the more realistic situation of 
finite size systems subjected to hydro dynamical expansion and consequent flow. 
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APPENDIX A: MASS-SHIFTED IN THE ENTIRE VOLUME 


It turns out that the formalism can be presented in the simplest manner if we assume, 
that the whole thermalized medium is filled with mass-shifted quanta, and that the whole 
medium decays suddenly to asymptotic quanta. This case is the subject of this appendix. It 
is also possible, that e.g. due to density inhomogeneities, the volume where the mass-shift 
is non-zero is different (smaller) than the totality of the volume filled out by thermalized 
quanta. This case will be investigated in the next appendix. 

For the case of non-relativistic hydrodynamics, assuming for the sake of simplicity a sud¬ 
den freeze-out, / = / d 3 r f dtS(t — t 0 )Eij = Eij f d 3 r, the chaotic and the squeezed 

amplitudes are easily obtained from the ones previously derived in Eq. © and © by 
taking the limit R s —► oo in those equations as well as in all the others that immediately 
follow them, i.e., Eq. ©, ©, and ©. The intuitive way to understand this limit is to 
consider the mass-shift region as extended so as to include the entire volume of the system, 
by simply taking the limit R s —> oo. The volume of the system, however, will be still delim¬ 
ited by the Gaussian profile with rms R. In this limit, the last two terms in Eq. © exactly 
cancel, since 

lim I{ 2 (Rs, R, (u), m) = /^(oo, R, (u), m). (Al) 

R s —>00 ’ ’ 

Consequently, the effective squeezing region becomes R s —> R. 

Before writing the resulting expression for G c (l,2), G s (l,2) and G c (i,i), it is uscfull to 
define two parameters in terms of which we can write those expressions, i.e., the flow-modified 
temperature T* and the flow-modified radius of the single volume case 


/?* = lim p(R s ,R,mu )) 

R s —>oo 


(A2) 


as given in Eq. ©• We can then write the chaotic amplitude, G c (ki, k 2 ) as 

( Ei 2 


G^ y (ki,k 2 ) = lim 


Tin 


R s —xx> ( (27 t) 3 1 ° 

s 0 | 2 (2vri? 2 ) 3 ^ 2 exp 


(|c 0 | 2 + |s 0 | 2 ) Iip(Rsi R, T, m*) + |s 0 | 2 /jj 2 (-R s , R, 00 , m*) j 




(2tt) 3 ' 

El, 2 * /| |2 1 I |2\ 

(J c ol + l s °l ) 


R 2 


(k]_ — k 2 )" 


x (27T.R 2 ) 3 / 2 exp 


(k 2 + k 2 ) Rl \ m(u)(ki + k 2 ) 

(ki — k 2 ) T i' 


4m* T 


n 2 


2 m*TR 


(A3) 
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We can rearrange the above terms in a more compact form, and explicitly writing in terms 
of the variables defined in Table HI (from which we can see that R 2 /RT = R/T.) we have 


Gl v ( k!,k 2 ) = ^l / 2 R 3 \s 0 \ 2 exp 


R 2 (lr , 

- — ( ki - k 2 


^(|c 0 | 2 + | S0 r) (2vri? 2 ) 3 ^ 2 exp 


( k l + k 2^ 
8 


x exp 


im(u)R f 2 


2 




(A4) 


TABLE I: Parameters used in the single volume case. 


Parameter 

Relation to other parameters 

Integral results where they appear 

T* 

T * = T + f>> 2 

I?j(oo,R,T, m.) 

R. 

R f- = R-i (i + agf) 

Ifj(oo,R,T, m.) 


The coherent amplitude G s ( k i 7 k 2 ) can be written as 


G^ v ( k i, k 2) = lim 


Ei,2 


R a ^oo ( ( 27 t ) 3 


Co |So| 


2 n* Q I{ 2 {Rsi R, T, m*) + I^ 2 (R S , R, oo, m*) 


El ’ 2 c 0 |s 0 |(27ri? 2 ) 3/2 


(2t r) 3 
x (2 tuR 2 ) 3/2 exp 


exp 


-^-( ki + k 2 ) 2 


+ 9ml s I 

+ 72 ^ 2co|So1 


(k? + k|) _g/ ( ki + k2 ) 2 

4m*T 2 l 2m.TR v 2) 


(AS) 


Similarly to what was done before, can also rewrite the expression for G s ( k i, k 2 ) explicitly 
in terms of the variables defined in Table IH leading to 

Et ' 2 : c 0 \s 0 \(2rrRV'\ 


Gy(k„k 2 ) = 


(2tt) j 


exp 


-^(ki + k 2 ) 


^ 3/2 ( 27r ^*) 3/2 (2 c 0 1s 0 1) exp 


( k i - k 2 f 
8 m.T 


x exp 


im(u)R 


( k i + k 2 ) ^ — 


8 ^ + fW^ 


(A6) 


2 m*T* 

Also, the single-particle distribution, the amplitude appearing in the denominator of both 
the BBC and the HBT correlation functions can be written as 




IV t 
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= lim 


E, 


l,l 


R a ^°O (27r) 3 L 

(^i s «i 2 < m2 ) 3/2 

Eijn, 


n o f|co| 2 + | Sol 2 ) Ii,i(Rs, R, T, m*) + \s 0 \ 2 Ili(Rs,R, oo,m*) 


(l-ol 2 + KI-) (2^P 


exp 


kf / m(u) k, t 

2m*T 2 l m*TR 


(AT) 


Analogously, we can rewrite G c (kj,kj) explicitly in terms of T* and R *, as 


^r(k.) = ^^(k^ko 

- (§ y ' s °' 2 < 2rf > 3/2 + 0 C »' 2 + ' s »' 2 ) < M -) 3/2 -- V 2m ,T.y 


' exp 


1 (A8) 


Let us investigate the vanishing squeezing and the large momentum limits of the single 
particle spectra, similarly to the analysis of the correlation functions as was done after Eq. 

P|). 

In case of vanishing squeezing, ?7i* —> m, |s 0 | —> 0 and |c 0 | —> 1, hence we spectra will 
contain a thermal and a flow contribution, and we recover the results of Ref. zj|. In the 
large momentum limit, for non-vanishing squeezing, rather surprisingly the single particle 
spectra becomes a constant. This corresponds to the decay of a modified vacuum with a fixed 
volume, described by the first term of (1A8I) . This is the direct consequence of our neglecting 
for the present purposes the position dependence of the in-medium mass modification. Also, 
this result implies that the squeezing mechanism not only makes strong signals in the back- 
to-back correlations, but there is also an interesting signal for squeezing in the single particle 
spectra. 


APPENDIX B: MASS-SHIFTED IN PARTIAL VOLUME 

If the mass-shift occurs only in a certain portion of volume V s of the whole system 
U(> V,), the expressions for the amplitudes contain other terms besides the ones discussed in 
the Appendix El Again, in order to avoid too much clutter it is useful to define appropriate 
flow-modified variables. However, in this case, we will need to define two sets of such 
parameters as flow-modified radii and temperatures, one set corresponding to the region 
where there is no mass-shift, which we will denote by TZ, p, and T, and another for the 
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inside of the mass-shifted region, denoted by R** and T*, this last one, naturally, being the 
same as defined in Tabled 


TABLE II: Parameters used in the two-volume case. 


Parameter 

Relation to other parameters 

Integral results where they appear 

T 

T = T + m(u) 2 

If j(oo,R,T,m) & I c (R s ,R,T,m ) 

n 

K~ 2 = R~ 2 (1 + %(u) 2 ) 

Ii,j(oo,R,T,m) 

p 

p~ 2 = 7Z~ 2 + R~ 2 

I i,j{ R s-,R, T ,™) 

T* 

T* = T+£(u) 2 

Iij(Rs,R,T, m*) 

P* 

p: 2 = n- 2 {i + ^) + Rf 

I?j(R s ,R,T, m*) 


In the case where the mass-shift occurs in a small portion of the system volume, V s < V, 
the chaotic amplitude is given by Eq. (ED and ED, be., 

-Ei ,2 


1,2) = j^k[n*o (M 2 + M 2 ) I{ i2 (RsiR,T,m*) + \s G \ 2 T l2 {R Sl oo, oo, m*) 
+ n 0 Ii 2 ( oo, -R, T, m) — n 0 I^ 2 (R s , R , T, m) 


(Bl) 


Working out each of the integrals above separately, we have 


I c h2 (R s ,R,T,m*) = (27rp2) 3/2 exp J ^ kl + 


4 m*T 


(ki - k 2 ) + i 


.m{u)( ki + k 2 ) 


2m-.TR 


= (27Tp2)3/2 exp 

/ 1 


,m(u)p2(k2-k2) 


—«- 


x exp 


_I P* l _ k. P_ 

8 m*T 2 y 1 2; 8m*(R 2 T + R 2 T*) 


2 m*RT 

2N J (-R 2 + -R 2 )(ki + k 2 ) 2 


/ p2 (-R s , oo, oo,m*) = (2tcR 2 s ) 3/2 


exp 


Ii 2 (oo, R,T,m) = (2nTZ 2 ) i R exp 

= (27r7t 2 ) 3/2 exp 

/ 1 




(k? + k 2 ) n 2 


(B2) 

(B3) 


4 mT 


(k]_ — k 2 ) + i 


,m(u)(k 1 + k 2 ) 


1 2 


2mTR 


\u)TZ 2 { k 2 -k 2 

* 2RT 


x exp 


— + 5b (kl _ k2> 2 _ 

8 mT 2 { ; 8 mT 


(B4) 


Ii 2 (R s , R, T, m) = (27rp 2 ) 3 / 2 exp 


M + k l 

4 mT 




(ki - k 2 )+i 


.m(u)( ki + k 2 ) 


2mT i? 
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= (27rp 2 ) 3/2 exp 


>)p 2 (k 2 -k 2 ) 
2 RT 


x exp 



(kx - k 2 ) 2 - 


(i? 2 + ^)(k 1 + k 2 ) 2 ' 
8m(R 2 T + R 2 f ) 


(B5) 


Substituting these into Eq. ED, the complete expression for the chaotic amplitude can 
finally be written as 


Gf(l,2) = n*pl\\c 0 \ 2 +\s 0 \ 2 

( 1 


x exp 


£|(k. 


exp 
2 \ 


8 m*T 


, m(u)p:{kj - k 2 ) 

2 m*RT 

(i? 2 + ^)(k 1 + k 2 ) 2 
2j 8 m,{R 2 T + R%) 
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!, 2 d 3 I „ 12 


(27t) 3 /2 


exp 


i? 2 


(kx - k 2 ) 


^{(2^)^exp 


— I 


. (w)7t 2 (k 2 — k 2 


2 i?T 


x exp 




\8 mT 


(k, - k 2 ) 2 - (kl + k2)i 


- (27rp 2 ) 3/2 exp 

( 1 


—*- 


»p 2 (k 2 -k 2 ) 


x exp 


_ r 

8 mT ^ 2 


2RT 


(kx - k 2 ) 2 - 


8 mT 


(R 2 + R 2 a )( kx + k 2 ) s 


8 m{R 2 T + R 2 T) 


(B6) 


The coherent amplitude is given by Eqs. © and © , as 
Gf (k t ,k2)= El ' 2 


2-c 0 |s 0 | 2-n.Q I[ 2 (R s , R, T, m*) + I^ 2 (R S , oo, oo, m*) , (B7) 

f/7TJ L J 


where 


Il 2 (R s ,R,T,m *) = (27rp 2 ) 3/2 exp 
= (2 t rp 2 ) 3/2 exp 
x exp 


(k? + k 2 ) p 2 


4m*T 


1 + i 


. m(w) 


2 mJTR 


(kx + k 2 


—i 


.m{u)pl( kx + k 2 ) 5 


2 m*RT 


(kx-k 2 ) 2 p 2 , , . 2 ( J R 2 + J R 2 )(k 1 + k 2 ) 2 


H. _—fki + kB 2 - 

8 m*T 2 K 1 V 8 m*(R 2 T + R 2 T*) 


I* 2 (-R s , oo, oo, m*) = (27ri? 2 ) 3//2 exp 


“Y (kx + k 2 ) 2 


, (B8) 


(B9) 


By substituting the above two terms, Eq. (El and (El into Eq. (ED, we get the final 
form of the coherent amplitude 


E\ 


Gl V {h 2 ) = jYyj (2ttR 2 s ) 3/2 (c 0 s 0 ) exp 


i? 2 


(kx + k 2 
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( 27r ^) 3/2 ( 2c o s o) exp 


—i 


• ™( M )P*( k i + k 2) 5 


x exp 


2 m*RT 
(R 1 2 + R 2 s )( kl + k 2 ) 2 


( k i ~ k ^) 2 _ d (k k)2 _ 

2 ' 1 2> 8m*(R 2 T + R?T< 


8 m*T 

The single-particle distribution for this situation of two volumes can be written as 


(BIO) 


JVfpk.) = Gf (i,i) 


2V( 


J l,l 


(2vr) 


n o (I Co | 2 + | So | 2 ) IZi(Ra,R,T,m *) + |s 0 | 2 /^(i? s , oo, oo, m*) 

( R 2 + W 


+ rc 0 /£(oo, f?, T, m) - n 0 Ii ti (R s , f?, T, m) 


^j|(2^)^| So | 2 + nj (2 x P D 3/2 (kol 2 + l»o| 2 ) exp 
n 0 (2Tril 2 ) 3 4 5 6 7 8 9 / 2 exp ^- 


, no(2vrp 2 ) 3/2 exp 

2 mT / 


2m*(R 2 T + R 2 Tt 
(R 2 + R s 2 ) k 2 


(BH) 


2m(R 2 T + R S 2 T) _ 

Note again that in the large momentum region only the first term survives, which cor¬ 
responds to a constant contribution given by the modified vacuum, and the size of this 
contribution is less than in case of Eq. (lA8j) . as the in-medium modified quanta do not fill 
the entire volume in the present case: R s < R. 
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